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Abstract 



We evaluate the two-loop corrections to Bhabha scattering from fermion loops in the context 
' of pure Quantum Electrodynamics. The differential cross section is expressed by a small number 

of Master Integrals with exact dependence on the fermion masses m e ,m/ and the Mandelstam 
invariants s,t,u. We determine the limit of fixed scattering angle and high energy, assuming 
the hierarchy of scales <C <C s,t, U. The numerical result is combined with the available 
non-fermionic contributions. As a by-product, we provide an independent check of the known 
electron-loop contributions. 



1 Introduction 



Bhabha scattering is one of the processes at e + e _ colliders with the highest experimental precision 
and represents an important monitoring process. A notable example is its expected role for the 
luminosity determination at the future International Linear Collider ILC by measuring small-angle 
Bhabha-scattering events at center-of-mass energies ranging from about 100 GeV (Giga-Z collider 
option) to several TeV. Moreover, the large-angle region is relevant at colliders operating at 1-10 
GeV. For some applications a full two-loop calculation of the QED contributions is mandator;^]. 

A large class of QED two-loop corrections was determined in the seminal work of [2]. Later, 
the complete two-loop corrections in the limit of zero electron mass were obtained in [3] thanks 
to the fundamental results of [4,5]. However, this result cannot be immediately applied, since the 
available Monte-Carlo programs (see e.g. [6-13]) employ a small, but non-vanishing electron mass. 
The a 2 ln(s/mg) terms due to double boxes were derived from [3] by the authors of [14], and the 
close-to-complete two-loop result in the ultra-relativistic limit was finally obtained in [15,16]. Note 
that the diagrams with fermion loops have not been covered by this approach. The virtual and 
real components involving electron loops could be added exactly in [17, 18]. The non-approximated 
analytical expressions for all two-loop corrections, except for double-box diagrams and for those with 
loops from heavier- fermion generations, can be found in [19] . For a comprehensive investigation of the 
full set of the massive two- loop QED corrections, including double-box diagrams, we refer to [20-22]. 
The evaluation of the contributions from massive non-planar double box diagrams remains open so 



In order to add another piece to the complete two-loop prediction for the Bhabha-scattering cross 
section in QED, we evaluate here the so-far lacking diagrams containing heavy- fermion loops. The 
cross section correction is expressed by a small number of scalar Master Integrals, where the exact 
dependence on the masses of the fermions and the Mandelstam variables s, t and u is retained. In 
a next step, we assume a hierarchy of scales, rr? e <C m?j <C s,t,u, where m e is the electron mass 
and m/ is the mass of a heavier fermion. We derive explicit results neglecting terms suppressed 
by positive powers of mg/m 2 , m^/x and TnrJx, where x = s,t,u. This high-energy approximation 
describes the influence of muons and t leptons and proves well-suited for practical applications. In 
addition, we provide an independent cross-check of the exact analytical results of [17] (we used the 
files provided at [23] for comparison) for mf = m e . 

The article is organized as follows. In Section [5] we introduce our notations and outline the 
calculation and in Section [3] we discuss the solution for each class of diagrams. In Section [4] we 
reproduce the complete result for the corrections from heavier fermions in analytic form and perform 
the numerical analysis. Section [5] contains the summary, and additional material on the Master 
Integrals is collected in the Appendix. 

2 Expansion of the Cross Section 

We consider the Bhabha-scattering process, 
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and introduce the Mandelstam invariants s, t and u, 
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1 Note that leading two-loop effects in the electroweak Standard Model were already incorporated in [1], 
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Figure 1: Classes of Bhabha-scattering one-loop diagrams. A thin fermion line represents an 
electron, a thick one can be any fermion. The full set of graphs can be obtained through proper 
permutations. We refer to [24] for the reproduction of the full set of graphs. 



where m e is the electron mass, E is the incoming-particle energy in the center-of-mass frame and 9 
is the scattering angle. In addition, s + t + u = 4 m 2 .. 

In the kinematical region m 2 <C s,t,u the leading-order (LO) differential cross section with 
respect to the solid angle reads as 
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?(t + t2 + st ) + h(i + s2 + st ) + h {s + t)2 



(2.5) 



where a is the fine-structure constant. At higher orders in perturbation theory we write an expansion 
in a, 

da _ da LO (a} dcr NLO fa\2dcr NNLi 
dn ~ ~ ''" ' 

Here dcr NLO and dcr NNLO 



d£l V 7r J dfl V 7r J dfl ^ ^ 

summarize the next-to-leading order (NLO) and next-to-next-to-leading 



order (NNLO) corrections to the differential cross section. In the following it will be understood 
that we consider only components generated by diagrams containing one or two fermion loops. 



2.1 NLO Differential Cross Section 

The NLO term follows from the interference of the one-loop vacuum-polarization diagrams of class 
la (see Figure [T]) with the tree- level amplitude, 
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Here (x) is the renormalized one-loop vacuum-polarization function and the sum over / runs 
over the massive fermions, e.g. the electron (/ = e), the muon (/ = /i), the r lepton (/ = r). Qf is 
the electric-charge quantum number, Qf = — 1 for lcptons. 

In this paper we will focus on asymptotic expansions in the high-energy limit. In order to fix our 
normalizations explicitly, we reproduce here the exact result for n^ 1 "* (x) in dimensional regulariza- 
tion. Adding n^ c *(a;), the counterterm contribution in the on-mass-shell scheme (see the following 
discussion in Subsection 12. 3p . to U^P un (x), the unrenormalized one-loop vacuum polarization func- 
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where e = (4 — D)/2 and I? is the number of space-time dimensions. The normalization factor is 



F r 



(2-11) 



/i is the 't Hooft mass unit and je is the Euler-Mascheroni constant. Standard one-loop integrals 
appearing in Eq. (|2.8p are defined by 
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Note that Master Integrals with 1 lines and an internal scale m were derived in [20,24] setting m = 1. 



For the present computation we introduce a scaling by a factor 21 and we get 



A (m f ) = F e — 5- Tlllm, 
\m 2 / 



B (x,m f ) = F e 



SE212m [x]. 



(2.14) 



(2.15) 



In the small-mass limit, Ao vanishes (the result for Tlllm can be read in Eq.(4) of [20]), and the 
one-loop self-energjQ reads as 
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+ 2L f (x) + -L 2 (x) 



(2.16) 



where we introduced the short-hand notation for logarithmic functions (in our conventions the 
logarithm has a cut along the negative real axis), 



L f (x) = In 



Finally, neglecting 0(m 2 /x) terms, Il^\x) reads as 
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Note that the O(e) term in Eq. (|2.18p is not required for the NLO computation, but it will become 
relevant at NNLO. Here n^ 1 '^) will be combined with infrared-divergent graphs showing single 
poles in the e plane for e = 0. The exact result for ILy(x) is available at [24]. 

2 Here, the argument x of SE212m [x] is one of the relativistic invariants s,t,u. This deviates from earlier conventions, 
where we denoted by x the dimcnsionless conformal transform of s, t, u. This remark applies also to Master Integrals 
in the Appendix. 
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Figure 2: Classes of Bhabha-scattering two-loop diagrams containing at least one fermion loop. We 
use the conventions of Figure [1] Note that class 2a contains three topologically different subclasses. 
We refer to [24] for the reproduction of the full set of graphs. 

2.2 Outline of the NNLO Computation 

At NNLO we have to consider: 

• The interference of the two-loop diagrams of classes 2a-2e (see Figure [5]) with the tree-level 
amplitude; 

• The interference of the one-loop vacuum-polarization diagrams of class la with the full set of 
graphs of classes la-lc (see Figured]). 

The complete result can be organized as 

j.nnlo J_2ixtree J^laxli 



2— loopxtree 1 — loopx 1 — loop 

In order to compute the NNLO differential cross section we use the following reduction strategy: 

• The generation of all the diagrams is simple and has been made with the computer-algebra 
systems GraphShot [25] and qgraf /DIANA [26-28]. We spin-sum the squared matrix elements 
and take the traces over Dirac indices in D dimensions using the computer- algebra system 
FORM [29] . The resulting expressions are combinations of algebraic coefficients depending on 
s, t, u, m e , mf and e and two-loop integrals with scalar products containing the loop momenta 
in the numerators. An example showing the complexity of the result (two-loop box diagram 
of class 2e, see Figure [2]) can be found at [24] . 

• We reduce the loop integrals to a set of Master Integrals by means of the IdSolver implemen- 
tation [30] of the Laporta algorithm [31,32]. The complete list of massive Bhabha-scattering 
Master Integrals can be found in [20] . 

Next, we evaluate the Master Integrals: 
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Figure 3: Counterterm-dependent Feynman rules relevant for Bhabha scattering for i — 1 (one loop) 
and i = 2 (two loops) . Note that in the on-mass-shell scheme e 2 = Ana at all orders in perturbation 
theory. 

• Integrals arising from graphs of classes la-lc (Figure[l}, 2a-2c (Figure^]) and 2d-2e (Figure^ 
with electron loops) have been computed exactly through the method of differential equations 
in the external kinematic variables and expressed through Harmonic Polylogarithms [33] or 
Generalized Harmonic Polylogarithms [34,35]. Here we agree perfectly with the work of [17, 
23] . Non-approximated results for the various components of the differential cross section are 
collected in a Mathematica [36] file at [24]. 

• Integrals generated by the diagrams of classes 2d-2e (Figure [2] with heavy- fermion loops) are 
computed through a method based on asymptotic expansions of Mellin-Barnes representations. 
We derived appropriate Mellin-Barnes representations [37, 38] for each Master Integral and 
performed an analytic continuation in e from a range where the integral is regular to the origin 
of the e plane [4,5]. This is done by an automatic procedure implemented in the package 
MB.m [39]. To proceed further, we assume a hierarchy of scales, m 2 e <C ml <C s,t,u, where 
/ =/= e. After identifying the leading contributions in the fermion masses (in the same spirit as 
in [40]), we express the integrals by series over residua, and the latter are sumed up analytically 
in terms of polylogs by means of the package XSUMMER [41]. Asymptotic expansions for the 
master integrals with two different masses were given in [42]. They, and also few lacking 
expansions of simpler masters needed here have been collected in Appendix [AJ We refer for a 
detailed discussion to [22], where the technique was employed to derive approximated results 
for the massive Bhabha-scattering planar box master integrals. All the mass-expanded masters 
may also be found in a Mathematica file at [24]. 

2.3 Renormalization 

In the following we will always deal with ultraviolet-renormalized quantities. After regularizing the 
theory using dimensional regularization [43,44], we perform renormalization in the on-mass-shell 
scheme. Here we relate all free parameters to physical observables: 

- The electric charge coincides with the value of the electromagnetic coupling, as measured in 
Thomson scattering, at all orders in perturbation theory; 

- The squared fermion masses are identified with the real parts of the poles of the Dyson- 
resummed propagators; 



5 



- Finally, field-renormalization constants are chosen in order to cancel external wave-function 
corrections. 



Counterterm-dependent Feynman rules are shown in Figure [3] Note that the presence of infrared 
divergencies at NNLO requires to compute one-loop counterterms including 0(e) terms. 



One-Loop Counterterms 

The one-loop counterterms read as 
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where the last equation follows from the U(l) QED Ward identity. In the ultrarelativistic limit, the 
one-loop fermion-mass counterterm is not needed, since it is always multiplied by the fermion mass. 
Note however that the same counterterm is relevant for the exact computation. 



Two-Loop Counterterms 

At the two-loop level we get 
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(2.23) 
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The result for SZ 2 ee is obtained including just fermion-loop diagrams and neglecting 0{m 2 /rri 2 ^) 



terms for / 7^ e. The expression for SZ 2 (as well as the one-loop counterterms of Eqs. (|2.20D - 



(|2.22p ). instead, is exact, since it follows from the single-scale diagrams of classes 2a-2b of Figure O 
Finally, we observe that the two-loop counterterm with two fermion lines is not required, since the 
use of an on-mass-shell renormalization removes external wave-function factors. 



3 Two-Loop Corrections 

In this Section we show our approximated results for all the components of the NNLO differen- 
tial cross section of Eq. (|2.6|) . Our short-hand notation for logarithmic functions can be found in 
Eq. (|2.17p . In addition, we define two combinations of the Mandelstam invariants: 

Vi(x,y;e) = x 2 + 2y 2 + 2xy ~ ex 2 , (3.1) 
V 2 (x,y;e) = (x + yf - e (x 2 + y 2 + xy), (3.2) 

where x(y) = s, t, u. Note that for e — these functions are proportional to the kincmatical factors 
appearing in the Born cross section of Eq. (|2.5[) and the NLO corrections of Eq. (|2.7[) . Moreover, we 
introduce a compact notation which will prove useful in discussing box corrections in Subsection [231 
and the complete NNLO differential cross section in Section 01 

L(Rf) = ^ H?J ■ (3-3) 
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3.1 Vacuum-Polarization Corrections 



The interference of the vacuum-polarization diagrams of classes 2a and 2b with the tree-level am- 
plitude can be written as 
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Here we introduced the auxiliary functions A (x) and A (x), which are expressed through the 



renormalized one- and two-loop vacuum-polarization functions IL^P(x) (see Eq. (|2 . 18[) ) and Uj^x), 



A 2a (x) = ^g)Re[n^ 2) (x) 
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where the result for 11^ (x) in the small fermion-mass limit reads as 
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Note that 0(e) terms in Eq 
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coming from the kinematical coefficients of Eq. (|3.ip can be safely 



neglected, since both I^p(x) and U^\x) are infrared-finite quantities. 



3.2 Vertex Corrections 

The contribution of reducible (irreducible) vertex corrections to the NNLO differential cross section 
can be readily derived from diagrams of classes 2c (2d) in Figure [H 
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Reducible diagrams 

The auxiliary functions A 2 ^ (x) and A 2 ^ (x) are given by the product of the renormalized one-loop 
vacuum-polarization function 11^ (x) (expanded in Eq. (|2. 18[) including 0(e) terms) and the renor- 



malized one-loop vector and magnetic vertex form factors Fy'(x) and i 7 ^ (x) 



Af{x) = ^2Q}Re[F^(x)Ilf\x) 
f 



I = V, M. 



The asymptotic expansion of Fy^ (x) is given by 



L e (x) 



1 + | C 2 - \ L e {x) - \ l\(x), 



(3.9) 



(3.10) 



whereas F^p(x) vanishes when we neglect the electron mass, F^'(x) — 0. The renormalized one- 
loop vertex develops an infrared divergency, which shows up as a single pole in the e plane for e = 0. 
Therefore, when computing the cross section, we sum over the spins the squared matrix element 
and we evaluate the traces over Dirac indices in D = 4 — 2 e dimensions. The needed kinematical 
structures include 0(e) terms (see Eq. (13. 
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Irreducible Diagrams 

The renormalized two-loop vertex diagrams of class 2d are free of infrared divergencies. Therefore, 
we can neglect 0(e) terms in the kinematical coefficients of Eq. (|3.1[) appearing in Eq. (|3.8[) , setting 



v a (x,y,;e) = v a (x,y;0), for a = 1,2. The auxiliary functions Ay d (x) and A^(x) contain the 
renormalized two-loop vector and magnetic vertex form factors (see [45-47] for a detailed discussion) , 



f 



I = V, M. 



(3.11) 



For the case with an electron loop, F Ie (x), the exact results in terms of Harmonic Poly logarithms, 
can be readily expanded in the high-energy limit. For the vector term we get 

?P)f^ _ 1 f 383 s.\ , 1 f 265 , s. \ r ^ , 19 r2 t ^ , 1 r 3, 



= 4(^7 - C 2 j + - 6 + C 2 j + -^(,) + (3.12) 

(2) 

For FyHx), / e, we perform an asymptotic expansion of the Master Integrals arising in the 
computation (see Table V in [20]) and we fully agree with the result of [48], 
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Since collinear logarithms are absent, the logarithmic structure of Eqs. (|3.12|) and (|3.13|) is obviously 
the same. 



3.3 Box Corrections 

The contribution of the renormalized two-loop box diagrams of class 2e is given by 
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Here the auxiliary functions can be conveniently expressed through three independent form factors 
B$ (x, y), where i = A, B, C, 
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Electron Loops 



(2) 

For the case with an electron loop, B\^(x,y), we get exact results in terms of Harmonic Polyloga- 
rithms and Generalized Harmonic Poly logarithms. An asymptotic expansion in the limit m 2 e -C s,t,u 
leads to 
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Heavy-Fermion Loops 

The list of Master Integrals here is given in Table V of [20]). At variance with the electron- loop case, 
it is not possible to compute them exactly by means of a basis containing Harmonic Polylogarithms 
and Generalized Harmonic Polylogarithms. Therefore, we use the high-energy asymptotic expansion 
discussed in Subsection 12.21 The results, expressed by the logarithms of the fermion masses L(Rf) 
(see Eq. (|3T3"j)). are: 



Bfj[x,y) 



1 2 
e 3 



2x 



- 5 
-3 



L(R f ) + L e {y) L e (x) 



9 



+ 



+ 
+ 



1*L 

3 y 



- ~L\R f ) + 



{2(— - IOC2 - 2C 3 ) " 2( y - 6C 2 )L(i2/) + gi 2 (i?/) 
82 



;) - 2Ca - 3L(i?/)ji e (x) - 2 



23 
~6~ 



I + 8C 2 - ii(i?/)]i e (y) 

5 .3, 



-2L(i2,) ^(y)+4 2-L{R f ) L e {x)L e {y)-A -L A e {y) 



- Le(x)L 2 M 



6 C2 + In" 



(I)M' 



12 



2 In 



L 3 (i?/) + 2 



121 



/ 262 

V~2T 
10 



9C 2 -4C 3 



() - yL(i?/)ji e (x) - 2 



- 4 

r 10 



(IW-I) 

- - 3C 2 )i(i2/) 



+ 12 C 2 - 2L(i2 / ) L e (y) 



13 



2L(fy) Lfc) 



r 16 



2L(fy) ^(2/) + 2 



17 



2L(fy) L e (x)L e (y) 



+ - Ll(x) + 6 L e (x) Ll(y) - 2 Ll(y)-2[6( 2 + In 2 In (l + |) 

- 4 K!)^H)^ Li 3H)} + i{ 2 (f-^- 2 cs) * 



- 2(^-3C 2 )i(fl/) 



V 9 



Q-L 2 (Rf) - -L 3 (Rf) 



6 + 12( 2 -3L(R f )\L e (y) + 
10 



LOR,) L^(x) - 



rl30 10 

1 ~v ~ y 

25 



£(i2/)j£e(s) 



+ 2[y - L{R f )\ L e (x)L e (y) + ^I%{x) - L'M + 3L e (x)L 2 e (y) 

6<2 + h2 (I)] h(l + |)_ 2ta (|) Ll2 (_|) +2Li3 (_|) } , 



(3.20) 



\\(^— +2x + y)\l- L ( ./? /■ ) h L,.(Y> 



e 3 v y 
2 x 2 r 262 
371 
112 



9 
23 



27 

2C2 



20 C2 
10 



3 

4Cs- ( 
L(R f )] L e (x) + 



50 

y 



12C 2 )i(i2/) 
2 



^(i?/) - - L 3 (Rf) 

16 C2 + i(i?/)] i e (y) 
r 5 



- 2L(R f ) Li(y) + 25- 2L(i2,) L e {x) L e (y) - 4 -^(y) 



- Mx) - [6C2 + In 2 (|) ] in (l + |) - 2 in (J) Li 2 (-f ) + 2Li 3 (-f ) } 

+ ^{f£ - 9C2 -4 C 3 - 2(| - 3( 2 )L(R f ) + 7 -L 2 (R f ) l -L 3 (R f ) 

Le(z) 



27 



9 

y — L(Rf) J -E^O^) 



9 

y + 12 C2 - 2L(i2 / )]L e (») 



20 i 

— - 2L(R f )\ L e (x)L e (y) 



+ -Ll(x) + 3L e (x)L 2 e (y) L 3 e {y) - [6(2 + 1* 2 (£)] In (l + f ) 



-(f) *H) 



2Cs 

1 

2 



/25 

ly 



3C 2 )L(i? / ) + ^i 2 (i? / )-^L 3 ( J R / ) 



6 



6 + 12C2 - 3L(i2/) L e (y) + 



l-L[R J )]L\{x)-\ 



65 5 

y ~ 3 

25 



L(R f )\L e (x) 



( . - )J -^e(y) 



10 



y - L(R f )]L e (x)L e (y) + ll*(x) - ^(y) + ^L e (x)L 2 M 
2 \x 



y 



In 1+i -In [Z. Li 2 



Li, - 



•!)}■ 



(3.21) 



^g/O^) = -7i77 i ~ L (%) + L -^) L e (x) + ^[x + y) ~ - L{R f ) + L e (y) L e {x) 



1 2 x 2 r 5 
e 3 y L 3 

| f! { _™ + 10<2 + 2c, + (| - . 6 )kjw - £j? (*,) + 

41 



(y - C2 - §£(*/)) M*) + (I + 8C2 - \L{Rf)) L e (y) 
- 2 (2 - L(R f ))L e (x) L e (y) + (§ - L(R f ))L 2 e (y) + ^L 3 e (y) - 2L e (x)L 2 e ( y ) 

+ N + 5 ■»■(!)]"»(' + 1) (-*)-"• (-!)}■ <" 2 > 

In order to study the numerical effects of massive leptons in two-loop box diagrams we consider the 
interference of the box diagram of class 2e (see Figure [2]) with the s-channel tree-level amplitude, 



B 



20,/ 



4 s 2 



: Re 



(3.23) 



where Baj can be found in Eq. (13.17j) for electron loops, and in Eq. (|3.20j) for / 7^ e loops. In 
Table Q] (Table [5]) we show numerical values for the finite part of B2 C j at values of y/s typical for 
meson factories, Giga-Z, ILC, and at two selected small and wide scattering angles, 9 — 3° (9 = 90°). 



B 2cJ [nb] / yfs [GeV] 


10 


91 


500 


e [see Eq. (|3.17|>] 


188758 


5200.08 


284.711 


fi [sec Eq. (|3.20|)j 


1635.62 


1686.88 


130.579 


T " 






39.5554 



Table 1: Numerical values for the finite part of i?2c./ of Eq. (|3.23p in nanobarns at a scattering angle 
= 3°. The first two entries for the r lepton are not shown since here the high-energy approximation 
in not justified (the same consideration applies to the top quark). 



B 2eJ [nb] / yfs [GeV] 


10 


91 


500 


e [see Eq. 1)3.170] 


143.162 


3.23102 


0.160582 


H [sec Eq. p.20pj 


61.3875 


1.79381 


0.0995184 


r " 


10.0105 


0.935319 


0.0639576 


t 






-0.00256757 



Table 2: Numerical values for the finite part of i?2e,/ of Eq. (|3.23|) in nanobarns at a scattering 
angle 9 = 90°. The first two entries for the top quark are not shown since here the high-energy 
approximation in not justified. 

For comparison we show in Figure the real part of the vertex function, see Eq. (|3 . 13[) . 
We see that the contributions from the box diagrams with heavier fermions are not strongly 
suppressed, but are instead of about the same size as the boxes with electron loop. This is different 
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[GcV] 


10 


91 


500 


e 


-124.237 


-254.293 


-400.574 




-4.8036 


-29.1057 


-70.1032 


r 




-2.08719 


-13.4901 



Table 3: The real part for the vertex form factor, see Eqs. (I3.12|) and (|3.13j) . 



to the self-energy and vertex corrections and may be traced back to the logarithmic structure of the 
contributions Eqs. (|3.20[) - (|3.22p . where terms of the order L^(x) appear. Further, in Eq. (|A.7[) we 
may see that this Master Integral has a dependence on L\ (x) , in contrast to the vertex and self-energy 
masters with heavy fermion loops. That originates in an additional collinear mass singularity from 
the external legs of this diagram diagram. One may control this easily by evaluating the singularity 
structure of the corresponding massless box diagram where only a scale M due to the internal loop 
exists, and see there some 1/e 2 terms which are absent in the corresponding SE and vertex diagrams. 
This leads finally to the fact that the two-loop corrections from heavier fermions are not numerically 
suppressed compared to the electron loop contributions. 



3.4 Products of One-Loop Corrections 

Finally, we consider the simpler components generated by the interference of one-loop diagrams 
among themselves. We start with the interference of diagrams of class la, 

J lax la 2 , i i 

= — {^i(M;0)A laxla ( S , S ) + -v 1 (t,s;0)A 1 ^(t,t) 

+—v 2 (s,t;0) U laxla ( S ,t) + A laxla (M)l )• (3.24) 
St L J J 

Here the auxiliary function A xla (x, y) contains the product of the renormalized one-loop vacuum- 
polarization function Ily (x) (see Eq. (|2.18[) ) with its complex conjugate, 



A laxla (x,y) = £ Q%Q%nf{x) [Uf 2 \y) 

h,h 

The interference of diagrams of class la with those of class lb gives 



(3.25) 



da 



lax lb 



a 2 r 1 
s I s z 



+ — 



v 1 (s,t- 1 e)A 1 v ^ lh (s, S ) + S 2 A^ lh ( S , S ) 
v 1 (t,s;e)A 1 v ay - lh (t,t) + 1*A%* lb (t,t) 



V2 



(s,t;e) (4 



axlb/ ,\ i a laX lb /± \ 

(s,t) + A w (t,s) 



+ l(s 2 A^ lh (s,t)+t 2 A^ b (t,s)) 



2 
2st 



( A la 



• lh {s,t) + < axlb (M))] }• 



(3.26) 



The auxiliary function A laxlb (x, y) is given by the product of Fy 1 ^ (x) and (x), the renormalized 
one-loop vector (see Eq. (|3.10j) ) and magnetic (vanishing in the high-energy limit) form factors for 
the QED vertex, and the complex-conjugate renormalized one-loop vacuum-polarization function 
n ( , 1} (x) (see Eq. (|2X8])1. 



A?* lh (x,y) = ^Rc{Fj (1) (x) [n^(»)]*}, I = V, M. 



(3.27) 
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Finally, the interference of diagrams of class la with those of class lc gives 



J laxlc J r i i 

-A\ Mc (s,t) + -Al^ lc (s,t) 

-S t 



4~7 



(3.28) 



Here the auxiliary functions A 1 axlc (s, t) and A\ ay ~ lc (s,t) take the form 

A\^ lc (s,t) =F e Y l Q 2 f^{[B^(s,t) + Bg\t,s) + 4 x) («,i) - B^( M )] [n^s)]*}, 
/ 



(3.29) 



Al^ lc (s,t) = F e ^(#Tte{[B%\a,t) + Bg\t,8)+B%\u,t)- Bg\u,a)] [< } ( S )]*}, 



(3.30) 



Al^ lc (s,t) = F e J2 Q}Re{[B^\ s ,t) + B<p(t,8) - B™(u,t) + B%\u,s)] [!#>(<)]*} 



(3.31) 



U^}\x) is given in Eq. (|2.18[) . and the new functions, in the small mass limit, read as 



B { l\x,y) 



-- (— + 2x + y) L e (x) + — I6C2 + 4L e (x) + 2L'i{y) 
e V y I y 



4L e (x)L e (y) 
2L e (x) L e (y) 
2L e (x) L e (y) 



2x 

y 



IOC2 + L e (x) + L e (y) - L 2 e (x) +Ll(y) 
IOC2+ 2L e (x) + 2L e (y) - L 2 e (x) + L 2 e (y) 



(3.32) 



--(2— + 2x + y )L e (x)+4— 8C2 + L'i(y) -2L e (x)L e (y) 
e V y J y 



2x 
II 



B ( c\x,y) 



IOC2- L e (x) + L e (y) -L 2 e (x) + I?M - 2L e (x)L e {y) 
IOC2+ 2L e (x) + 2L e (y) - L\(x) + L 2 e (y) - 2 L e (x) L e (y) 

Ax 2 - 2 r 



-L e (x) + 2- 

e y y 

4(.t + y)L E {x). 



8C2- 2L e {x) - L 2 e {y) + 2L e {x)L e {y) 



(3.33) 



(3.34) 



For the computation of the non-fermionic corrections these functions are needed up to first order 
in e, since they are combined with the real emission. However, this higher-order expansion is not 
relevant here. 



4 The Net Fermionic NNLO Differential Cross Section 

In this Section we use the results of Section [3] and derive an explicit expression for the NNLO 
differential cross section of Eq. (|2.19|) . 
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Note that the full set of two-loop fermionic virtual corrections to Bhabha scattering represents 
an infrared-divergent quantity. In order to obtain a finite quantity, we take into account the real 
emission of soft photons^) from the external legs of one- loop fermionic diagrams (class la, Figure [1]). 
The exact result is available in the literature, see e.g. Eq. (25) and Appendix A in [18]. Here we 
show the high-energy approximation relevant for our computation. We consider events involving a 
single soft photon carrying energy to in the final state, 



(pi) 



(p 2 ) 



(pa) + e+ (pi) +7(*), 



(4.1) 



and compute one-loop purely-fermionic corrections. Obviously, these real corrections factorize and 
their structure is completely equivalent to the tree-level ones. In complete analogy with Eq. 
we write 



where 



da 1 



2 da^° 



(4.2) 



dtt 

da^° 

dn 



2 s 



2 vi(s,t',e) 



2^"i(M;e) + j^v 2 {s,t;e) F(u>, s,t,m 2 e ), 



(4.3) 



(t,s;e) Q) Refn^(t) 



-l- t v 2 {s,t;e) Q/Re[n^( s ) + nf'ft)] } F{u, a, t, m 2 e ). 



(4.4) 



II^^x) can be read in Eq. (|2.18|) and, at variance with Eqs. (|2.5|) - (12.7I) . the kinematical factors 
introduced in Eq. (|3.ip need to be expanded up to 0(e), since the real-emission factor shows an 
infrared divergency, 



F(ui, s, t, ml 



2 
e 

In 2 
4 In 

4(2 



In 



m 
75 



2Li 2 - 



In 



2 In 



In [ — ] - In [ 1 



In 2 1 



2Li 2 1 



-In [ I 
2 In 



1 



2ui 

t 
s 



In 



- 1 



In 1 



(4.5) 



Summing the virtual contributions of Eq. (|2.19D to the real-photon emission of Eq. (|4.4p we write 
the NNLO fermionic corrections to Bhabha scattering through the sum of electron-loop contributions 
(dtr NNLO ' c ) and components arising from heavier fermion loops, 



do 



NXLO d(T NLO 



do~^ 



dfl 



dn 



</<> 



Jr-NNLO.f 



dn 



Of- 



da 



NNLO.f 



f*e 



dn 



J^NNLO^f 

E Q%Ql^^-(^ 

fuh^e 



dn 



3 The energy u) carried by a soft photon in the final state is small with respect to the center-of-mass energy E 
introduced in Eq. 12. 2t . 
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The double summation over the fermion species arises from the loop-by-loop terms of Eqs. (|3 . 6|) 
and (|3.24p . Here we do not include the case f\ = f% = e, which is incorporated in <i<7 NNLO ' . Note 
also the term proportional to Q^, coming from Eq. (|3.5[) . The result for electron loops can be found 
in Eq. (46) of [18]. For heavier fermion loops we introduce x = —t/s and get: 



da 



NNLO.f 4 



d(7 NNLO,2f 



dn 



da 



NNLO.f 



dn 



NNLO.f 



+ 
+ 



+ 



a 2 ( 
2s I 



{l-x + x 2 Y 



In (£} + ln( % ) + 4< 3 - i] + Hz)(± -S + l-i)}- < 47) 



+ x 2 )' 



3x 2 

Hx)\\n(R fl 



ln ( — ) + HRh)HRh) + In ( — ) (ln(fl A 



In 2 (a;) 



ln(i?/ 2 )-y 



21n( ^ 

mi 



4 7 a; 

3a; 6 3 
/ 1 11a; 
V3x2 ~ 2a; + 2 ~ 6 



)}■ 



NNLO.f 



« r n 



(l - a: + a; 2 ) 



NNLO.f i 

a, 'In 



(*)}■ 



In [ -4 



3a; 2 L 3 

589 37 



hr [ ^ 



+ HRh)-f) 

5 + 4a; - 2x 2 ) 

(4.8) 
(4.9) 

ln (1 - x) - ln {x) 



mi 
4795 



— + y In (R f ) - In 2 (%) - 2 ln (R f ) (ln (x) - In (1 - a;)) - 8Li 2 (a 



409 

~ls~ 



108 

40 t- / n! , ( s 
y L l2 (x)} + ln - 



19 



ln (R f ) + y In 2 (R f ) - In 2 (R f ) (in (x) - ln (1 - a:)) - 8 ln (R f ) Li 2 (a;) 

In 2 m (__L + il 

v ' V 3x 2 12a; 



2 4 11 23 16 , 

C2 k H 1 x H a; 2 

1 3a; 2 3a; 2 3 3 



5 

4 ~ 
4 

3x 



x 

12 
1 

" 2 



2 2 A , 2 ,., s / 2 11 5 
+ r ) +ln (1 ^H-3^ + fe-2 + 



In(x)( 



55 
Vite 2 



11 

6a; 
83 

9x 



5 
2 

65 
IT 



11 

Y a 

85 
18 : 



-*(* 
V3a; 2 

10 

3^ 



31 iQ 

6a; 

rln 3 (l- 



- -x 2 ) +ln(x)ln(l - x) 
o / 

+ y^ 2 ) +|ln(l-x) (- 
31 10 ,M 1 , , , , / 1 31 11 x x 2 \ 

x 2 \ , , , , , s ( 1 1 4 



1 

3^ 



1 

.7; 



4 
3 



x 2 \ 1 , , , , 9 , > / 1 2 7 a;\ , o , , 
y ) + -ln(, ) ln 2 (l-, ) (-- + -- i + -) + ln 2 (, ) 



55 



46 



ln(.Rf) ^ H h -a;' 

1 /; V 3a; 2 12a; 4 12 3 

, , n s ! 2 11 5 11 2 2 



20 2 , /„ x S 2 4 1 5 
+ Y X +ln( ^H3^-3x-2 + 3 



a; a; 

3 



+ ln 2 (l-a;) 
+ ln (a;) ln (1 — x) 
+ In (a;) 



56 



27 



19 



-.x 2 +Li 2 (.x)(-- + --7 + 3a;--a; 2 )+hi(i? / )( 7 ^-- + ---a;+- 



14 

18X 2 ~ 9x + Y 
10 29 9 29 

9a- 2 ~ Ete + 2 ~ ~9 X 
10 37 1 

9x 2 + l8^ + 2 
589 1753 701 

54x 2 + 108a; ~ ~36 
37 56 47 67 



x 

4 10 

— X — ^ 

9 9 



10 

9 

25 



— x 



9 
925 

To8 : 

10 



X 

10 



10 



2x' 



31 10 31 



+ ln(l -x) 
10 



56 
27a; 2 



9a; 2 
161 



9.r: 
56 



6 
161 



18 



9 



56 
— i 
27 



\n(R f ) x H f x x 

1 j> 1 9a; 2 18a; 3 18 9 



54x 9 54 
, 2 20 32 20 , 
+ C2 ( -— + — - — + — x - 2x 



X" 



3.c 



Li 3 (a) ( 



V3x 2 



7 

3.7' 



3-^a; 
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+ (2 



19 13 43 311 98 , , ,„ . / 2 4 11 23 16 ,\ 



3 

x 



C3(-3^ + r-5 + -»-2^ 



11 

y 



(4.10) 



NNLO.f 



(l - x + x 2 )' 



In 



+ ln(%) -ln(l -a;) + In (x) In (J?/) 



+ 
+ 4 



--ln(%)] [l + m(l-z)]} 
" / 4 7 



In 



■ I ln ( x ) [ Tl. ~ T~ + 6 ~ q x + ~i x ' 
mi I \6x z 6x 3 3 



3 — —x 



( 2 1 
ln{x) ln(l-x)^ - - 



1 1 



i 2 

\^-x 
16 29 „ 23 



In 2 (a;) 



1 1 

- + 1--X 



3 lnW b-^ +13 -y s+ r)l- (4 - n 



In order to have compact results we used 

(— l)" +p_l 



S n , P (y) 



(n-l)!p! J 



1 , ln"- 1 (x)ln P (l-^) 
ax ■ — . 



(4.12) 



In Table 2] (Table [5]) we show numerical values for the NNLO corrections to the differential cross 
section for a scattering angle 9 = 3° (9 = 90°). In both tables we set u> = E/10. Finally, in Figure |4] 
we plot the ratio of the two-loop fermionic corrections to the tree-level cross section, 



R(VJ> - (H) 

for V« = 10 GcV and ^ = 500 GeV. 



2 acr NNLO + da^° 
da LO 



(4.13) 



da 1 dQ [nb] | v 


s [GeV] 


10 


91 


500 


LO QED 


[Eq. 


440873 


5323.91 


176.349 


LO Zf itter 


[49,50] 


440875 


5331.5 


176.283 


NNLO (e) 


[Eq. gl])] 


-1397.35 


-35.8374 


-1.88151 


NNLO (e + fi) 


-1394.74 


-43.1888 


-2.41643 


NNLO (e + /x -1 


■r) " 






-2.55179 


NNLO photonic 


[14,16] 


9564.09 


251.661 


12.7943 



Table 4: Numerical values for the NNLO corrections to the differential cross section respect to the 
solid angle. Results are expressed in nanobarns for a scattering angle 9 = 3°. Empty entries are 
related to cases where the high-energy approximation cannot be applied. 

It is clear from the Tables, that although there is no decoupling of the heavier fermions (as 
indeed there shouldn't, since the typical scale of the process is large compared to all the masses), 
the electron loop contributions dominate in the fermionic part and the latter is still substantially 
smaller than the pure photonic corrections. 



5 Summary 

In this article, we completed the computation of the virtual two-loop QED fermionic corrections 
to Bhabha scattering. Based on the kinematics of the targeted phenomenological applications, we 
considered the limit m\ <C ml <C s, t, u. 
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do- / dfl [nb] | ^/s [GeV] 


10 


91 


500 


LO QED [Eq. 


0.466409 


0.00563228 


0.000186564 


LO Zf itter [49,50] 


0.468499 


0.127292 


0.0000854731 


NNLO (e) [Eq. JH])] 


-0.00453987 


-0.0000919387 


-4.28105 • 10~ 6 


NNLO (e + fi) 


-0.00570942 


-0.000122796 


-5.90469 • 10- 6 


NNLO {e + n + t) " 


-0.00586082 


-0.000135449 


-6.7059 • 10~ 6 


NNLO {e + fi + t+ t) " 






-6.6927 • 10" 6 


NNLO photonic [14, 16] 


0.0358755 


0.000655126 


0.0000284063 



Table 5: Numerical values for the NNLO corrections to the differential cross section respect to the 
solid angle. Results are expressed in nanobarns for a scattering angle 9 = 90°. Empty entries are 
related to cases where the high-energy approximation cannot be applied. 

| , 1 , 1 , 1 n | , 1 , 1 , 1 n 



-0.025 



-0.01 



-0.02 



-0.03 




50 



100 



150 



Figure 4: Ratio of the fermionic NNLO corrections to the differential cross section respect to the 
tree-level result for y/s = 10 GeV and s/s = 500 GeV. A solid line represents the electron-loop 
contributions, a dotted one the sum of electron- and muon-loop ones, and a dashed one includes also 
t leptons. 



The fermionic double box contributions with two different mass scales have been derived for the 
first time here. Their numerical importance is comparable to the two-loop self-energies and vertices. 
We note, however, a qualitative difference. Due to the structure of the collinear singularities of the 
graphs, the contributions of the heavier fermions are not suppressed. 

A numerical estimation of differential cross sections shows that the net fermionic two-loop effects 
may be neglected for applications at LEP 1 and LEP 2, but have to be taken into account for 
precision calculations when a level of 10 -4 has to be reached, as is anticipated for the Giga-Z option 
of the ILC project. 

Completing the NNLO program for Bhabha scattering requires still several ingredients. First, let 
us mention the contributions from the five light quark flavors. Here, an approach based on dispersion 
relations a la [51] should be suitable. On the other hand, the heavy top quark might be considered 
decoupling in a large part of the interesting kincmatical regions. Furthermore, an implementation 
of the loop-by-loop corrections with pentagon diagrams has to be done. 

Finally, light fermionic pair emission diagrams need to be considered. As known from the form- 
factor case, they are responsible for the cancellation of the leading part of the logarithmic sensitivity 
on the masses. 

Exact and approximated results are made publicly available at [24]. The combination of our re- 
sult with the photonic two-loop corrections of [16] and with electron loop corrections of [17, 23] 
proves well-suited for phenomenogical purposes, e.g. a precise luminosity determination at a future 
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Figure 5: Same as Figure [H including the photonic contributions of [2, 14, 16] (dash-dotted lines). 
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Note added 

We would like to thank T. Becher and K. Melnikov for drawing our attention to a problem with a first 
version of our result, which lead us to discover an incorrectly expanded integral. After correction, 
Eq. 4.9 agrees with the result published in the meantime in [53]. 



A Mass-Expanded Master Integrals 



The list of Master Integrals required by our computation can already be found in Table V of [20] . 
The eight most difficult masters, those involving two different mass scales, have been derived in [42]. 
Because they are a substantial part of the present study we reproduce them here: 



SE312Mlm[on shell] = 



SE312Mlmd[on shell] = 



M 2 m- 4£ 



R 



If 



- - -L(R) 
18 3 v ; 



1 16 4^ 3 4 v ; 



iL 3 (i?))+i? 2 (-J + ^(i?)-iL 2 (i?))]}, 



-4e 



I 2e 2 + 2e 



2L{R)\+-{l+C,2) + L{R) 



L\R) 



L 2 (R) 



(A.l) 
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V412Mlm[x] = 



V412Mlmd [x] 



V412M2m [x] 
V412M2md [x] 
B512M2m [x,y] 



B512M2md [x,y] 



+ 



+ 



~ (3 + 3C 2 - 2C 3 ) + (1 + Ca) L (R) + L 2 (R) + \h 3 (R) 
o 3 



R 



1 



L(R) 



l-L(R) + h 2 (R) 



R z 



■S + 5 £ W + <-S + S^ + I L '^)]}' 

19 - 3C 2 - L 2 Jx) 



(A.2) 



1 



1 



-4e f 

™ lie" ' 2e ' 2 

Af 2 r 

2 + 4C 2 - 4C 3 - 2L m (x) + 2L M {x) - 4C 2 Lm(x) 

x L 

+ 2L m (x)L M {x) - L 2 M (x) - L m (x)L M {x) + L 3 M (x) }, 



(A.3) 



m r 1 1 

to 2 1 2e 2 e 
M 2 r 1 1 

- 1 + 3(2 + im(^) + L M (x) 



1 + + 2 - C2 + L m {x) + - L 2 m (x) 



L m (x)L M {x) 



-4ef 1 

m < — k 
l2e 2 



-L 



M 



(*)]}, 



(A.4) 



G.r 



+ L m (xj\ + ^(19 + C2) + 5 L m {x) + L 2 n (x)}, (A.5) 

(A.6) 



12Cs - 6(2Lm(x) - L 3 M (x) 

1 . , s 1 



= ^— |^L m (.x) + -(-C 2 + 2L m (x) + -L 2 n (x) + L m (x)L m (y)^ 

- 2( 2 - 2( 3 + AL m {x) + L 2 m {x) + h 3 m (x) - 4( 2 L m (y) 
+ 2L m {x)L m {y) + L m {x)L 2 n {y) - ^L 3 n (y) 

- (3C2 + ~L 2 m (x) - L m (x)L m (y) + h 2 m (y)) In (l + |) 



fi TO (a:) - L m (y]jLi 2 



(A.7) 



xy 



4e r 1 r 1 

— <^ L m (x)L m (y) + L m (x)L(R) - 2( 3 + ( 2 L m (x) + AC, 2 L m (y) 

11 I e I J 



2L m {x)L 2 m (y) + -L 3 n (y) - 2( 2 L(R) + 2L m (x)L m (y)L(R) - -L 3 {R) 
(3C2 + \lUx) - L m (x)L m (y) + \L 2 m {y)) In (l + f ) 



(A.8) 

We list also the other expanded masters, including the correct normalizations. Note that, compared 
to the conventions employed in [20] and in Eq. (|2.16[) , all integrals are rescaled by a factor m L ( D ~ 2l \ 
where L is the number of loops, D = 4 — 2e and I is the number of internal lines. Expansions are 
performed up to the order required by our computation. For example, we include 0{m 2 ) terms in 
SE212m[x] (see Eq. (fATTO]) ) since the reduction procedure generates coefficients containing 1/m . 
The same consideration applies to 0(e) terms, which are included as long as the reduction brings 
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inverse powers of e in the coefficient functions. Since in the following no ambiguities arise, in we 
drop the subscript / and we set L(x) = h\{—m 2 /x), 



Tlllm = (m z ) 



2\l-e 



l + f (l + |)+ e 2 fl+ C2 



(A.9) 



f 1 777 ' 

SE212m[x] = mT 2t |- + 2 + L(x) + 2 — |^1 - L{x) 

+ Il 2 (x) + ^(2 + 2C 2 -L 2 (x))] + 

+ 4L(x)-l-( 2 L(x) + L 2 (x) + ]:L 3 (x) + —(2 + 0. 

+ 4Cs + C2 L(x) - |i 3 (x))]}, 



4-f + 2L(x) 



-C2- 3C3 



(A.fO) 



SE210m[x] = m- 2e ^ + 2 + L(x) + e 4 - ^ + 2 L(.x) + 1 L 2 (x) ] 

+ £ 2 [8-C2 - ^Cs + 41(a) - ±C 2 £(aO + L 2 (x) + 1 L 3 (x) ] }, (A.ll) 



V3 llm[x] = !^-!l{4C 2 + ii 2 (x) - e[-5C 3 + C 2 i(x)-ii 3 (x)]}, (A.12) 



/■ 2x1 2e\ 1 5 ff 5 / 55 25 11 \ 

SESllm [on shell] = (m 2 ) 1 ^— + _ + _ + - C 2 + e (-- + y C2 + y Cs ) 



2 / 949 55 , 55 „ 303 A \ 



(A.13) 



SE312m[x] = ("> 3 ) 1 " 2e {^ + ^(3-^) + 5 + ( 2 -L 2 (x)- 



13 1 r / \ 

IT + 2^ 



+ e 



16 



3 + 3C 2 + ^Cs - 4L(x) + 2C 2 i(x) - 3L 2 (x) - L 3 (x) 



a; /115 (2 13 
m 2 V~L6~ ~ T + T 



L(x) + \l 2 (x))]}, 



(A.14) 



SE312md[x] = m -*{^ + ±-l + ^-L(x)-±L 2 (x)-^[-2 + L 2 (xj 

+ *[~ + ^ 2 + ^ 3 - + C2L ^) - 2L2 ^) - \ L *W 

2 

- ^(-6Cs - 4L(x) -2C 2 £(x) + + i 3 W)]}- 



(A.15) 



Finally, the mass expanded one- loop box Master Integral B412m [x , y] can be collected from Eqs. (4.70)- 
(4.75) of [52]: 



B412m[x,y] = 



rn 



' 2e (2 r 

-\-\m 

u lei 



xy 1 e 



-In I - 



+ 2L 2 (y)-2L(y)ln( - 

\y 



4C 3 - 9C 2 i(y) 
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+ -L 3 (y) + 5( 2 ln 



-L 2 (y)ln 



2 In - 



+ 2 In 



In 

Li 2 



In 1 



x 

y 

2LL, 



x 

y 



ln d 



In 1 



- 6C2 In 1 



x 



(A.16) 
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